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Abstract 

In the present work the linearized problem of plasma wave reflection from a 
boundary of a half-space is solved analytically. Specular accommodative conditions 
of plasma wave reflection from plasma boundary are taken into consideration. Wave 
reflectance is found as function of the given parameters of the problem, and its 
dependence on the normal electron momentum accommodation coefficient is shown 
by the authors. The case of resonance when the frequency of self-consistent electric 
field oscillations is close to the proper (Langmuir) plasma oscillations frequency, 
namely, the case of long wave limit is analyzed in the present paper. Refs. 17. Figs. 
6. 
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1. BASIC EQUATIONS 

Research of degenerate electron plasma behaviour, processes which take 
place in plasma under electric field, plasma waves becomes more and more 
actual at the present time in connection with the problems of such intensively 
developing microelectronics and nanotechnologies PQ [6]. 

This work continues the research of the electron plasma behaviour in 
external longitudinal alternating electric field [7] - jllj . 

In the present work linearized problem of plasma wave reflection from 
a boundary of a half-space of conductive medium is solved analytically. 
Specular accommodative conditions of electron reflection from plasma boun- 
dary are taken into consideration. The diffuse boundary conditions were con- 
siderated in [B]-[TU]. 
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Expression for wave reflectance is obtained and it is shown that in the 
case when normal electron momentum accommodation coefficient takes on a 
value of zero the wave reflectance is expressed by the formula obtained earlier 
in 0, HQ. 

Let us consider degenerate plasma which is situated in a half-space x > 0. 
We assume that self-consistent electric field E(r,t) inside plasma has one 
x-component and varies along the axis x only: E = {E x (x, t), 0, 0}. 

In this case the electric field is perpendicular to the plasma boundary 
which is situated in the plane x = 0. 

47re 2 iV 

Here oj„ is Langmuir (proper) plasma oscillation frequency, uj v = , 

m 

N is the electron numerical density (concentration), rn is the mass of the 
electron. 

Let us take the system of equations which describes plasma behaviour. As 
the kinetic equation we take r-model Vlasov — Boltzmann kinetic equation: 

at Or op r 

Here / = f(r,v,t) is the electron distribution function, e is the electron 
charge, p = mv is the electron momentum, m is the electron mass, r is the 
characteristic time period between two collisions, f eq = f eq (r,t) is the local 
equilibrium distribution function of Fermi and Dirac, f eq = 0(£p(£, x) — £), 
where Q(x) is the function of Heaviside, 

B(x) = 

8,p(t, x) = ^mvp(t, x) is the disturbed kinetic energy of Fermi, £ = ^mv 2 is 
the kinetic energy of the eletron. 

Let us consider the Maxwell equation for the electric field 

divE(r,£) = 47re J(f(r, v, t) - /o(v)) cM Fl . (1.2) 

where 

d^F = d3p = dpxd Py dpz - 

Here fo is the undisturbed electron distribution function of Fermi and 
Dirac, /o(S) = 0(£f — £), ^ is the Planck constant, £i? = \mvp is the 
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undisturbed kinetic energy of Fermi, vp is the electron velocity on the Fermi 
surface which is considered as spherical. 

Let us search for the solution of the system (1.1) and (1.2) in the following 
form 

v 

f = /o(£) + £ F 5(£ F - £)H(x, /i,t), /i = — . 

v 

We obtain (see P) linearized system of equations of Vlasov — Maxwell: 

l 

- + Mt^ + H(x 1: /i, ti) = fie(xi,ti) + 2 J H ( x ii tf-> *i)^ , > ( 1 - 3 ) 



dt 



de(xi,ti) 3cj^ 2 '' 



dx\ 2v 



2 



H(xi, //, tijdjj! . (1-4) 



-l 



Here e(x\,t\) is dimensionless function 

e{x,t) = ——E x {x,t), 

Vtf 

x\ = x/l is the dimensionless coordinate, where / = v f t is the average 
free path of electrons, t\ = vt is the dimensionless time , v is the effective 
frequency of electron scattering, v = 1/r. 



Supposing that k is the dimensional wave number, we introduce the dimen- 
sionless wave number k\ = k— : then we have kx = 1 1 , where e = — . 
Let us introduce the quantity loi = lot = to /v. 

2. BOUNDARY CONDITIONS STATEMENT 

Let the plasma wave move to the plasma boundary situated in the plane 
X\ = 0. The electric field of the wave changes according to the following law 

e+(xi, ti) = E x exp(-z( h uJih)). (2.1) 

£ 

The amplitude of this wave E\ we assume to be given. On the plasma 
boundary this wave reflects and the electric field of the reflected wave has 
the following form 

e-(x h ti) = E 2 exp(i( (2.2) 

£ 
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The amplitude E2 is unknown and is to be found from the problem 
solution. The quantities uj\ and k\ are not independent, the following depen- 
dence uj\ = uj\{k\) is determined from the solution of the dispersion equation 
which will be introduced below. 

It is required to determine what part of the wave energy is absorbed under 
the wave reflection from the plasma boundary, and what part of the energy 
is reflected, and also to find the phase shift of the wave. It means we have 
to calculate the reflectance which is determined as square of module of the 
ratio of reflected and incoming waves amplitudes 



R(k, cj, e) = 



E2 u 



(2.3) 



and to find the argument of the amplitudes ratio 

/E \ 

(f)(k, uj, e) = arg = arg E 2 - arg E x . (2.4) 

Let us outline the time variable of the functions H(xi, /i, t\) and e(xi, ti), 
assuming 

H(xi,nM) = e-^/i^i,//), e{x u t 1 ) = e^^e^i). (2.5) 

The system of equations (1.3) and (1.4) in this case will be transformed 
to the following form: 

1 

dh If 

H— + (1 - iui)h(x h 11) = fie(xi) + -J h(x h //)d//, (2.6) 

-1 

d ^§jH^W. (2.7) 

-1 

Further instead of X\, t\ we write x, t. We rewrite the system of equations 
(2.6) and (2.7) in the form: 

1 

dh If,, 

fj,— — h zoh(x, fi) = fie(x) + - / h(x, fi')dfi\ zq = 1 — iu\. (2.8) 



-1 



1 



d(ix) 3 ///(.r. //V//. (2.9) 



dx 2e 2 J 

-1 
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We consider the external electric field outside the plasma limit is absent. 
This means that for the field inside plasma on the plasma boundary the 
following condition is satisfied: 

e(0) = 0. (2.10) 

The non-flowing condition for the particle (electric current) flow through 

the plasma boundary means that 

1 

J fjLh{0,fj,)dfjL = 0. (2.11) 

-i 

In the kinetic theory for the description of the surface properties the 
accommodation coefficients are used often. Tangential momentum and energy 
accommodation coefficients are the most-used. For the problem considered 
the normal electron momentum accommodation under the scattering on the 
surface has the most important significance. 

The normal momentum accommodation coefficient is defined by the follo- 
wing relation: 

a p = JTZJTi <Ka p <l, (2.12) 
where P{ and P r are the flows of normal to the surface momentum of incoming 

to the boundary and reflected from it electrons, 

o 1 

P l = J n 2 h{S),n)dii, P r = J id 2 h(0,id)did, (2.13) 
-1 o 
quantity P s is the normal momentum flow for electrons reflected from the 

surface which are in thermodynamic equilibrium with the wall, 
l 

P s = J ii 2 h s (v)d^ r fl e h s (fi) = A 8 , < fi < 1. (2.14) 
o 

The function h s (ji) is the equilibrium distribution function of the correspon- 
ding electrons. This function is to satisfy the condition similar to the non- 
flowing condition: 

o 1 
J nh{Q,n)dn + J fjih 8 (fjt)dfjL = 0. (2.15) 

-1 o 
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We are going to consider the relation between the normal momentum 
accommodation coefficient a p and the diffuseness coefficient q for the case of 
specular and diffuse boundary conditions which are written in the following 
form: 

h(0, fi) = (1 - q)h(0, —fi) + a s , < // < 1. 

Here q is the diffuseness coefficient (0 ^ q ^ 1), a s is the quantity 
determined from the non-flowing condition. 
From the non-flowing condition we derive 

o 

a s = — 2q J fih(0, fi)dfi = qA s . 
-l 

After that we find the difference between the flows 

o 1 

Pi — P r = q J /i 2 /i(0, fi)dfi — J fi 2 a s dfi = 
-1 o 
o 1 
= Q J l^ 2 h(0, fi)dfi — q J /i 2 A s d/i = qPi — qP s . 
-1 o 
Substituting the expressions obtained to the definition of the normal momen- 
tum accommodation coefficient we obtain that a p = q. 

Thus, for the specular and diffuse boundary conditions the normal momen- 
tum accommodation coefficient a p coincides with the diffusion coefficient q. 

Together with the specular and diffuse boundary conditions other variants 
of boundary conditions are used in the kinetic theory. 

In particular, accommodative boundary conditions are widely used. They 
are divided into two forms: diffuse accommodative and specular accommo- 
dative boundary conditions (see |12]). 

We consider specular accommodative boundary conditions. For the func- 
tion h these conditions will be written in the following form: 

h(0, n) = h{0, -n) + A 1 + A 2 /i, 0<//<l. (2.16) 

Coefficients A\ and A^ can be derived from the non-flowing condition and 
the definition of the normal electron momentum accommodation coefficient. 
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The problem statement is completed. Now the problem consists in finding 
of such solution of the system of equations (2.8) and (2.9), which saatisfies 
the boundary conditions (2.10)— (2.16). Further, with use of the amplitudes 
of reflected and incoming waves found it is required to find the reflectance 
of the incoming wave energy (2.3) and the argument of the amplitudes ratio 
(2.4). 



3. THE RELATION BETWEEN FLOWS AND BOUNDARY 

CONDITIONS 

First of all let us find expression which relates the constants Ao, A\ from 
the boundary condition (2.13). 

To carry this out we will use the condition of non-flowing (2.12) of the 
particle flow through the plasma boundary, which we will write as a sum of 
two flows: 

1 o 
Nq = J /ih(0, fi)d/i + J /ih(0, n)dji = 0. 

o -1 

After evident substitution of the variable in the second integral we obtain: 

l 



h{0,fjL)-h{0,-fjL) 



d\i = 0. 



N = J fi 

o 

Taking into account the relation (2.16), we obtain that Aq = —2Ai/3. 
With the help of this relation we can rewrite the condition (2.16) in the 
following form: 

h(0, n) = h(0, -n) + A^ji - < ii < 1. (3.1) 

o 

We consider the momentum flow of the electrons which are moving to the 
boundary. According to (3.1) we have: 

Pi = Pr~ (3.2) 

It is easy to see further that 



P s = A s /3. 



(3.3) 
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With the help of the formula (3.2) we will rewrite the definition of the 
accommodation coefficient (2.12) in the form: 

u P Pr ~ « P y + ^(1 - a p ) = 0. (3.4) 

Let us consider the condition (2.15). From this condition we obtain that 

o 1 
A s = -2 J /ih(0, \i)d\i = 2 J fih(0, -/j,)d/i. 
-1 o 
Using the condition (3.1), we then get 

l 

A s = 2 J nh(Q, rfdiJL. (3.5) 
o 

Now with the help of the second equality from (2.13) and (3.4) we rewrite 
the relation (3.3) in the integral form: 

l 

(V - -^h(0,fi)dfi = -gg(l - (*p)M- (3.6) 

o 

Now the boundary problem consists of the equations (2.8) and (2.9) and 
boundary conditions (2.10), (3.1) and (3.6). 



4. SEPARATION OF VARIABLES AND 
CHARACTERISTIC AL SYSTEM 

Application of the general Fourier method of the separation of variables 
in several steps results in the following substitution: 

hr,(x, ii) = exp(-— )$(t/, fi), e v (x) = exp(-— )E{r)), (4.1) 

rj r\ 

where 77 is the spectrum parameter or the parameter of separation, which is 

complex in general. 

We substitute the equalities (4.1) into the equations (2.8) and (2.9). We 

obtain the following characteristic system of equations: 

1 

z (ri - fj)^(ri, fj) = W E(r]) + $(ry, fjTjdfjf, (4.2) 

-1 
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z E(rj) = -^-^ J <K^W- (4-3) 

-l 

Let us introduce the designations: 

UJ 2 £ 2 Z .1 + 7 2 
7= 1, m = — -, 2: = 1-2 , = 277^0. 

cj p 3 £ 

Substituting the integral from the equation (4.3) into (4.2), we come to 
the following system of equations: 

1 

(77 - f*)*(ri, fi) = ^( W - rft\ -rftE{rj) = \ [ $(ry, //) dp!. (4.4) 

zo * J 

-1 

Solution of the system (4.4) depends essentially on the condition if the 
spectrum parameter 77 belongs to the interval — 1 < 77 < 1. In connection 
with this the interval — 1 < 77 < 1 we will call as continuous spectrum of the 
characteristic system. 

Let the parameter 77 6 (—1,1). Then from the equations (4.4) in the 
class of general functions we will find eigenfunction corresponding to the 
continuous spectrum: 

= F(r,,»)^, (4.5) 

zo 

where 

F(??5 M) = p^L^l _ C M S ^ - M ). (4.6) 
77-/7 77 

In the equation (4.6) ^(x) is the delta-function of Dirac, the symbol Px~ l 
means the principal value of the integral under integrating of the expression 
the fuction X(z) is called as dispersion function of the problem, 

1 

X(z) = 1 + - I ^^d/i. (4.7) 
cj /I- z 

-1 

The function (4.5) is called eigenfunction of the continuous spectrum, 
since the spectrum parameter 77 fills out the continuum (—1, +1) compactly. 
The eigensolutions of the given problem can be found from the equalities 
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(4.7). The dispersion function X(z) we express in the terms of the Case 
dispersion function: 

X(z) = l-- + -(l-^)x c (z) : 

where 



1 1 

dr 1 f rdr 



T — Z 2 J T — Z 

-I -1 



is the Case dispersion function 

The boundary values of the dispersion function from above and below the 
contour are calculated according to the Sokhotsky formulas 



A ± (^) = A(/i) ± ^-{r,l - /i 2 ), — 1 < // < 1, 



from where we have 

Z7T 

A+(//) + A-( M ) 



A + (» - A {/I) = -g— /i(r/ 2 - ^i 2 ), 



where 



= A(^), — 1 < /i < 1, 



l 

A(^) = l + ^-/ ! ^^, 

2?7i^o J »7 — A* 

-l 



and the integral in this equality is understood as singular in terms of the 
principal value by Cauchy. Besides that, the function X(fi) can be represented 
in the following form: 

AW = l-- + -(l-4) *cto> A c (/,) = 1 + £ In 



5. EIGENFUNCTIONS OF THE DISCRETE SPECTRUM AND 

PLASMA WAVES 
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According to the definition, the discrete spectrum of the characteristic 
equation is a set of zeroes of the dispersion equation 

\(z)/z = 0. (5.1) 

We start to search zeroes of this equation. Let us expand take Laurent 
series of the dispersion function: 

X(z) = X 00 + ^l + - A +--- , |s| >1. (5.1) 

Z L 



Here 



Aoo EE A(oo) = 1 - — + 

A 2 = -I(i 



1 1 27 + is + 7(7 + is) 



z 3z r]j (1+7 + ^e) 2 

1/1 1 \ 9 + 5ie(l + 7 + ie) 



3b\3 brjfJ 15(1 + 7 + ^) 2 

1/1 1 \ 15 + 7ie(l + 7 + ie) 



A4 z U 777^ 35(1 + 7 + ie) 2 

One can easily see that in collisional plasma (i.e. when e > 0) the coefficient 
Aoo 7^ 0. Consequently, the dispersion equation has infinity as a zero r)i = 
00, to which the discrete eigensolutions of the given system correspond: 
hoo(x 7 n) = /i/zo, eoo(x) = 1. 

This solution is naturally called as mode of Drude. It describes the volume 
conductivity of metal, considered by Drude (see, for example, [14]). 

Let us consider the question of the plasma mode existance in details. We 
find finite complex zeroes of the dispersion function. We use the principle of 
argument. We take the contour r+ = Tr U j £ , (see Fig. 1), which is passed 
in the positive direction and which bounds the biconnected domain Dr. This 
contour consists of the circumference {Tr : \z\ = R, R = l/e,e > 0}, and 
the contour j £ , which includes the cut [—1, +1], and stands at the distance 
of e from it. 

According to the principle of argument the number [T5] of zeroes in the 
area D £ equals to: 

N = — idlnX(z). 
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Fig. 1. 



Considering the limit in this equality when e — > and taking into account 
that the dispersion function is analytic in the neighbourhood of the infinity 
we obtain that 

2m J d A-(r) " m J d A"(r) " n &lg A"(r) ' (5J) 

-10 



Consider the curve 

7={^: z = G(r), < r < +l}, 
where G(r) = A + (r)/A _ (r). It is obvious that G(0) = 1, limG(r) = 1. 

T— >1 

Hence, according to (5.3), the number of zeroes N equals to the double 
number of turns of the curve 7 round the point of origin, i.e. N = 2x(G), 
h{G) = Ind [0>+ i]G(r). 
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Fig. 2. 

Thus, the number of zeroes of the dispersion function belonging to the 
complex plane out of the segment [—1,1] of the real axis equals to the double 
index of the function G(r), calculated on the semi-axis [0, +1]. 

We take on the plane (7, e) (see Fig. 2) the curve L, determined by 
equations: 7 = — 1 + y / Li(r), e = y / L2(r), ^ r ^ 1, where 

T M 3^[ S 2 + Ao(1 + Aq)] 2 3tV 
llTj A [ S 2 + (1 + A ) 2 ] ' 2[T) Ao[ S 2 + (l + A ) 2 ]- 

Here 

A oW = l + Il„i^, s(r) = \r. 

In the same way, as was shown in the work [TB], we can prove that if 
(7, e) G D + , then the index of the problem equals to unity, i.e. N = 2 is the 
number of zeroes which equals to two, and if (7, £) £ D 7 then the index of 
the problem equals to zero, i.e. N = 0. 

We would like to notice, that in the work |T7] (see also [TB]) the method 
of examination of the boundary mode in the case when (7, e) G L was 
developed. 

Since the dispersion function is even its zeroes differ from each other by 
sign. We designate these zeroes as following ±770, by t]q we take the zero which 
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satisfies the condition Re 770 > 0. The following solution corresponds to the 
zero 770: 

ZqX , 



h m (x,/j) = exp 



Z X^E 2 T] fl 



m 



e m {x) = exp 



)E> 



2- 



(5.2) 



r/o z 770 - fi 770 

This solution is naturally called as mode of Debay (this is plasma mode). 
In the case of low frequencies it describes well-known screening of Debay 
[3]. The external field penetrates into plasma on the depth of t\d, Td is the 
raduis of Debay. When the external field frequancies are close to Langmuir 
frequencies, the mode of Debay describes plasma oscillations (see, for instance, 

mm)- 

From the equalities (2.2) for the wave e-(x,t) with the help of (2.6) and 
the equality (4.1) follows the relation between the wave number k and the 

zero of the dispersion function 77o(cj, v): i — = , hence 770 = 770(7, e) = 

1 + 7 .£ 

— + 1 F 

The equalities (2.5) and (2.6) jointly with (5.2) mean that the reflected 
wave corresponds to the zero 770: 

(kx 



rr _ Vl ~ W 

.ojiq — . . exp 



V £ 



cut 



e m = exp 



kx 

ut 

£ 



and the wave incoming to the plasma boundary corresponds to the symmetric 
zero —770: 



n -Vo / 1 \ ex P 



(kx 
i ^ h ut 



-1 10 



= exp 



(kx N 
% ^ h ut 



6. EXPANSION IN THE TERMS OF EIGENFUNCTIONS 

In the work [TT] it was shown that from the non-flowing condition (2.11) 
and the condition on the electric field (2.10) it results that the trivial (equal 
to zero) solution of the present problem corresponds to the point r]i = 00. 

We will show that the system of equations (2.8) and (2.9) with the boundary 
conditions (3.1), (3.5) and (2.10) has the solution which can be represented 
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as an expansion by the eigenfunctions of the characteristic system: 

E 2 riofjL-ril ,.kx E 1 r] ii + r\\ .kx 

h{x,fi) = ^-expu — ) H -exp(— % — )+ 

z 7/0 - fi £ z r/o + Li ' £ 

1 

H — / exp ( - z -)F(r],n)E(r])dr], (6.1) 
o 

i 

e(x) = E2 exp(i — ) + E'iexp(— i — )+ / exp ( — zq— )E(ri) dij. (6.2) 

£ £ J \ TjJ 



Here E\ is given, and E2 is unknown coefficient. Both of the variables 
(amplitudes of Debay) correspond to the discrete spectrum, E{rj) is unknown 
function, which is called eigenfunction of continuous spectrum. 

Our purpose is to find the coefficient of the continuous spectrum and the 
relation which connects the coefficients of the discrete spectrum. 

Let us substitute the expansion (6.1) into the boundary condition (3.1). 
We get the following equation in the interval < fi < 1: 

1 

1 2 
F(rj,fi)-F(ri,-ii)\E(ri)dri+(Ei+E2)Mv) = + (6.3) 

b 

Here 

<PoW = + : • 

id - r/o v + m 

Extending the function E{rj) into the interval (—1,0) evenly we transform 
the equation (6.3) to the following form: 

1 

2 

F(rj, fi)E(rj)drj + (E 1 + E 2 )<po(fi) - z Am = --zqAi sign /i. (6.4) 

-1 

Let us substitute the eigenfunctions of the continuous spectrum into the 
equation (6.4). We obtain singular integral equation with Cauchy kernel in 
the interval (—1,1): 

1 

(E 1 + E 2 )M») + [ V -^—^E(ri)dri - c^-E(jj) - z A lf i = 

J rj-n Li 

-1 



/ 
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2 

= --z Ai sign /i. (6.5) 



7. SOLUTION OF THE SINGULAR EQUATION 

We introduce the auxiliary function 

l 

^"^ EMd-n, (7.1) 
rj - z 

-l 

the boundary values of which on the real axis above and below it are related 
by the Sokhotsky formulas: 

M + {fi) - M-(fi) = 27n{fi 2 - rjl)E(fi). (7.2) 

M + (u) + M-(ll) , N 
^ ^ = M( M ), -1 < M < +1, (7.3) 

where 

l 

M(fj) = [ ^—^E(rj)dri : 

J 1 ~ V 

-l 

and the singular integral in this equality is understood as singular in the 
sense of principal value of Cauchy. 

With the help of the Sokhotsky formulas for the dispersion and auxiliary 
function we reduce the equation (6.5) to the boundary condition of the 
problem of determination of analytic function by its jump on the contour: 

^ + (aO + (Ei + E 2 )(fo(n) — z Aifj] — A~ (//) [M~ (/i) + (Ei + £ , 2 )^o(m) ~~ 

-zoAifi] = —^A 1 fi(r]j - fi) sign //, -1 < /i < 1. 

This equation has general solution (see [15]): 

l 

A(z) (/?(2;)(Ei + ^2) + M(z)-2;oAi2; = — — / d/z + Ciz, 

6c J fi — z 



-i 



where C\ is an arbitrary constant. 
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Let us introduce auxiliary function 



= 1 I M(M 2 -^) S ignM rf 

Then from the general solution we can easy find M(z): 

M(z) = -(Bi + S 2 )^(z) + z Aiz + gZo^i^y + -^4 



(7.4) 



Let us eliminate the pole of the solution (7.4) in the infinity. We get that 

C\ = — zqAiX^. 

Poles in the points z = ±r/o can be eliminated with the help of one equality 
since the functions constituting the general solution are uneven: 

(E 1 + E 2 )X'(r ]0 M-r ] i) 



ZqAi = 



(7.5) 



(2/3)T(t?o) - Xoolo 
We substitute the expansion (6.1) for the function h(x,/j) to the integral 
boundary condition (3.5). We get the following equation: 

l 



£?im(-7/o) + E 2 m(r] ) + J m(r))E{r))dr) = -^z A^ 



In (7.6) the following designations were introduced: 



1 — a r 



(7.6) 



1 1 

m(±77o) = J (V 2 - ^) F(±7]o, /i)dn, 771(77) = J (ji 2 - F(r), n)d/i. 

o o 
The coefficient of the continuous spectrum we will find from the Sokhotsky 
formula (7.1) after the substitution of the general solution (7.4) into it: 



E{rj) = 



1 



27n(?7 2 — rj\) 



2 fT+(ri) T~(ri) 
3 



1 



1 



A+(r7) A-(t7) 



(7.7) 



A+(r7) A-(t7) 

Let us notice that under the transition through the positive part of the 

cut (0, 1) functions T(z) and X(z) make jumps, which differ only by sign. 

Indeed, let us represent the formula for T(z) in the following form: 

l r -, 

1 1 

dji. 

Ii — z ii -t- z 

o 



i 

T(z) = Z -j(S- vl) 
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This integral can be calculated easily in explicit form. On the cut this 
integral is calculated according to the following formula : 

T(rj) = ^ [l + {rj 2 - rj 2 ) In - l)] , -1< r] < +1. 

Now from the Sokhotsky formula for the difference of boundary values we 
obtain that under condition < rj < 1 the following equality takes place: 

i7TT)(r) 2 — rj 1 ) 



\ + (r l )-\-(r ] ) = \(r ] )± 

T + (r ] )-T-(r ] )=T(r ] )± 
Now one can naturally find that 



c 

,2 Jl^ 



mr\{r\ — rj{) 



T + (r,)\-(r,) - T-(rt)X + (rt) = 2(T( V ) + Afa)) 



inr){rf — rj 2 ) 



With the help of the last relations we find the coefficient of the continuous 
spectrum from (7.5): 



E(rj) = z AiQ(r)), where Q{rj) = 3 



We introduce the integral 



Uz) = - f r J-iA 

c J rj - z 
o 



h[T(rj) + A(ry)] - A^^ 



(7.8) 



dr\. 



It is evident that in the complex plane this integral is calculated by the 
following formula: 



T (z) = * c [l + z + (z>- V l)\n(l-l) 



With the help of this function we represent the dispersion function in the 
form: X(z) = 1 — zTq(z) + zTq(— z), the function T(z) we also express in terms 
of this integral: T(z) = zTq(z) + zTq(-z). The sum of two last expressions 
equals to: \{z)+T(z) = 1 j t2zTq(—z). Let us note that the integral T(—z) is 
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not singular on the cut < 77 < 1. The sum X(rf) + T(r}) on the cut < r\ < 1 
is calculated in explicit form without applying to integrals: 

A(t?) + T(r]) = 1 + - \r) - 2t? 2 + 2r](r] 2 - rj{) \n(l/r) + 1) 

c L 

We calculate the integrals m(±?7o) and 771(77) m explicit form. The integrals 
m(±77o) can be determined easily: 

r 1 2 1 

m(±77 ) = -(r)l - r)l) - - + 770 + 770(770 - -) ln( 1) 

L 6 3 770 

Let us find the integral 777.(77). We have: 
1 



-^-277^0^-^(77)^(77). 

— 77 6 



Here #+(77) is the characteristic function of the interval < 77 < 1, i.e. 

« + w = { 1, 0<,?<1 ' 

L 0, -1<77<0. 

Computating the integral in the preceding equality we obtain that the 

integral 771(77) is calculated by the formula: 

1 2 r 

ra(»7) = (g - ^O? 2 - ^7i) + (^7 - 3) [ - c + 2t/ 2 - 7/(t7 2 - rff)f + (rj) 



where 



I + 77 

In -, < 77 < 1, 



In 



1 77 
1 — 77 

77 



-1 < 77 < 0. 



Now the equation (7.6) with the help (7.8) we rewrite in the form: 



Eim(-rio) + E 2 m(r]o) = -z Ai 



1 ~ Qfr 
36a„ 



m / 5 



(7.9) 



where 



Qm = J 777(77)Q(77)d77. 


After that substituting the variable Z0A1 into the equation (7.9) according 
to (7.5), we derive the equation: 



Eim{-j] {) ) + E 2 m(rio) = 



1 — a 



36a, 



P - + Q 



E\ + E 2 



m / 2 



T{t]q) - AooT/Q 
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from which we can find the amplitude required E^. 

£ 2 = a P m{-r]o)A{r)o) + g(?yo)C(^) 

We introduced the following designations in the formula (7.10): 

2 1 — a 

A(f] ) = -T(rj ) - AooT/o, C(a p ) = ^ p + a p Q m , 

B(r ]0 ) = (r ] l-r ] 2 )X , ( m ). 

Thus, all the coefficients of the expansions (6.1) and (6.2) are determined 
unambiguously, and this completes the proof of these expansions. 

It is seen from the equality (7.10) that under condition a p = we have: 
E2 = —Ei, i.e. under condition of pure specular reflection of electrons from 
the boundary the wave reflectance equals to unity: R = 1, and the phase 
shift of the incoming and reflected waves is equal to 180°, i.e. (j) = n, from 
where we see that arg E2 = arg E\ + 7r. 

Let us represent the formula (7. 10) in the form which is more convenient 
for numerical analysis. Let us designate the ratio of amplitudes as K, K = 
E 2 /E h the 

K= lH a P [m(r]o) - m(-r]o)] 

a P m(r]o) + C(a p )D(7/ )' 



where 



A(r]o) (2/3)T(t7o) -Aoor/o' 



8. LONG WAVE LIMIT 



For study of the incoming wave reflectance R = \K\ 2 and the phase shift 
(f) = arg K we will use the formula (7.11). 

We consider the dispersion equation with small values of the wave number 

We assume the frequency to is complex: u = luq + zcji. Then the quantity 
7 is complex also: 7 = 70 + *7i- Here 70 = uJq/uj p — 1, 71 = uj\/uj p . From 



21 



the equation (8.1) we find that when k is small 70 = 0.3k 2 , 71 = — 0.5e. We 
express the parameters of the problem in terms of k and e: Xoo = 0.6k 2 (l—i£), 

1 1 + 0.3/c 2 2 e, „ oj2n 1 + 0.3/c 2 + i0.5e 

z o = 2 " — ' 7/1 = -^(i + O- 3 ^ )> ^0 = 7 • 

With the help of these parameters let us carry out the study of the 
reflectance and the phase shift in long wave limit (when k is small by magni- 
tude). 

On the Fig. 3 one can see the dependence of the reflectance R on the 
wave number k for the case when e = 1CT 2 . The curves 1, 2, 3 correspond to 
the following values of the accommodation coefficient a p = 0.1, 0.5, 1.0. The 
curve 4 corresponds to the diffuse boundary conditions (see |H3]). From the 
graph it is seen that when k takes on small values the curve 3 (corresponding 
to dp = 1) coincides practically with the curve 4 (to which corresponds 
the case q = 1), which was obtained by the linearization of the reflectance 
value by k. Taking into account that under a p = the reflectance is equal 
to reflectance for specular boundary conditions (i.e. when q = 0) we can 
conclude that specular accommodative boundary conditions approximate 
specular and diffuse boundary conditions under a p = q very well. 

On the Fig. 4 the dependence of the reflectance R on the quantity e for the 
case a p = 1 is represented. The curves 1, 2, 3, 4, 5 correspond to the following 
values of the wave number k = 0.001, 0.01, 0.05, 0.1, 0.2. The more accurate 
analysis shows that with the growth of the accommodation coefficient the 
value of the reflectance decreases. 



On the Fig. 5 the dependence of the reflectance R on the value of the 
accommodation coefficient a p for the case e = 10 -3 is presented. The curves 
1,2,3,4 correspond to the following values of the wave number k = 0.05, 
0.1, 0.15,0.25. 

On the Fig. 6 the dependence of the angle (j) = oxgK (of the phase 
shift) on the quantity e for the case k = 0.2 is represented. The curves 
1,2,3 correspond to the following values of the accommodation coefficient 
a p = 0.1, 0.5, 1. The analysis shows that the dependence between the values 
of the angle (f) and the wave number and the accommodation coefficient as 
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Fig. 5. Fig. 6. 

well is small. 

9. CONCLUSION 

In the present work new boundary conditions for the questions of plasma 
wave reflection from the plane boundary of a half-space of degenerate plasma 
were proposed. These boundary conditions are naturally called as specular 
accommodative conditions. Such boundary conditions are most adequate for 
the problems of normal propagation of plasma waves (perpendicular to the 
boundary), since accommodation coefficient under such boundary conditions 
is normal electron momentum accommodation coefficient. 
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In the present paper the analytical solution of the problem of plasma wave 
reflection from a boundary with normal electron momentum accommodation 
is obtained. The analysis of the main parameters of the problem in long 
wave limit is carried out. This analysis shows that the boundary conditions 
proposed are intermediate between pure specular and pure diffuse boundary 
conditions. Indeed, from the Fig. 3 it is seen that all the graphs showing 
the dependence between the reflectance and the wave number are located 
between graphs corresponding to specular and diffuse boundary conditions. 
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